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Theoretical approaches to the photoionization of few-electron atoms are discussed. These include nonequi- 
Ubrium Greens functions and wave function based approaches. In particular, the Multiconfiguration Time- 
Dependent Hartree-Fock method is discussed and applied to a model one-dimensional atom with four electrons. 
We compute ground state energies and the time-dependent photoionization by the field a strong laser pulse with 
two different frequencies in the ultraviolet. 
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1 Introduction 

In recent years ultrashort, high-harmonic generated vacuum and extreme ultraviolet (vuv/xuv) laser pulses have 
become available. They open the way towards time-resolved observation of electronic dynamics in plasmas, 
atoms and condensed matter, for an overview see e.g. [1]. Electronic motion in atoms and relaxation processes 
can now be studied experimentally on sub-femtosecond time scales, e.g. [2^|, for a recent overview see |5|. This 
creates the need for complementary theoretical investigations. Interaction of laser radiation with matter has been 
actively studied in recent years using a quantum kinetic theory which was derived using nonequilibrium Greens 
functions (NEGF) methods, e.g. [6-8]. In contrast, application of NEGF to single atoms has only recently been 
attempted [9, 10]. While the ground state properties were very well reproduced, time-dependent calculations 
describing the ionization dynamics of atoms during and after a laser pulse turned out to be computationally very 
expensive [10]. 

In this paper we give a brief overview on the application of NEGF to time-resolved laser-atom interaction and 
outline the current problems. Furthermore, we discuss an alternative approach which is based on the A^— electron 
wave function - the Multiconfiguration Time-Dependent Hartree-Fock method (MCTDHF). The latter is quite 
efficient in application to few-electron atoms and small molecules. This is demonstrated by numerical results for 
a four-electron one-dimensional model atom subject to a short laser pulse. 

2 Theoretical concepts 

The unperturbed hamiltonian of N electrons in the atom is given by (we use atomic units) 

i=l 

where the potential of the nucleus, u(rj), as well as the two-particle Coulomb interaction w(rj—rj) = |ri — rj|~^ 
are assumed to be spin-independent. For times t > 0, the atom is disturbed by a time-dependent external field. 
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and the hamiltonian is modified: 



N 



H{t) = Ho- E(t)^r,, 



(2) 



j=i 



where the electromagnetic field is described in dipole approximation, and the field envelope is assumed to have a 
Gaussian shape with a pulse duration r and <o ^ ^ 



E{t) = 8ocos[uj{t — to)] exp 



{t - to) 
2t2 



21 



(3) 



2. 1 Quantum kinetic theory for a statistical ensemble of atoms 

To treat the interaction of an ensemble of atoms with the electromagnetic field a statistical theory is required. It 
can be based on density operators [11] or nonequilibrium Greens functions. Here we discuss the latter concept. 
Starting point are the fermionic field operators of the electrons (the nuclei are assumed to be fixed) und t/i, 
which guarantee that the spin statistics theorem is obeyed, 

V'(1)V(2) + ^(2)V(1) = V'^(l)^^(2) + V't(2)V't(i) = o, 
V^(1)V;^(2)+V^(2)V'(1) = <5(l-2), 

where ti = t2 and we denoted 1 = (ri, ii, Szi) [spin indices will be suppressed in the following]. The nonequi- 
librium properties of the electrons are completely determined by the two-time correlation functions (averaging is 
over the A^-particle density operator of the system) 

5>(l,l') = -Wl)V'ni')), 5<(ia')---(V'^(l'Ml)), (4) 

The two independent functions g> and g"^ contain the complete statistical information, in particular the nonequi- 
librium one-body density matrix, according to 

p(ri,r'i,i) = -*5<(l,l')lu=ti. (5) 

The equations of motion of and g< are the Keldysh/Kadanoff-Baym equations (KBE), e.g. [12, 13] which are 
the most general starting point for developing a quantum statistical theory of matter in a strong electromagnetic 
field. In recent years the KBE have been successfully solved numerically for spatially homogeneous plasmas, 
e.g. [14, 15] and for a weakly inhomogeneous electron plasma [16]. 

Yet these solutions require to propagate the correlation functions in the two-time plain which is computation- 
ally very expensive. Therefore, approximate methods have been developed. One possibility consists in expressing 
the correlation functions by the simpler single-particle distribution function f{t) [or density matrix p{t), Eq. (jsjl] 
via the so-called generalized Kadanoff-Baym ansatz [17]. For plasmas interacting with a strong field this ansatz 
has been generalized in refs. [6, 7]. The emerging quantum kinetic equation for the Wigner distribution which 
includes nonlinear field effects, inverse bremsstrahlung, collisional absorption and harmonic generation has been 
solved in ref. [8]. The second direction of recent research in the kinetic theory of laser-plasma interation was de- 
voted to include neutral atoms and their excitation and ionization. Starting from the KBE atomic Bloch equations 
were derived [18, 19] which allow to compute the absorption spectrum of an ensemble of atoms in thermodynamic 
equilibrium. This approach is very useful for weak excitation near the linear response regime. 

However, in the case of strong field ionization of atoms we expect essential nonequilibrium dynamics effects 
related to ionization of several electrons and electron-electron correlations developing rapdily on femtosecond 
and even sub-femtosecond time scales. These phenomena which have already been observed in experiments 
require a fully nonequilibrium approach. Starting again from the KBE we have in recent years developed a 
quantum kinetic approach for the electron dynamics in a single atom. The main difference to the previous 
applications to plasmas is that now the system of electrons under study is spatially strongly inhomogeneous - 
due to the confinement by the potential v of the atomic nucleus. To solve this problem the two-time correlation 
functions (|4]) were expanded into a basis of single-particle states or Hartree-Fock orbitals. This approach was 
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successfully tested for "artificial atoms", i.e. electrons in quantum dots in refs. [20]. More recently it was applied 
by the present authors to electrons in small atoms [10] such as a one-dimensional beryllium model atom following 
earlier work by van Leeuwen and co-workers [9] . In these works it was found that the ground state of few-electron 
atoms is well reproduced with simple approximations for the selfenergy such as the second Born approximation. 
The situation is more complex for time-dependent calculations. While lowest order calculations - using a time- 
dependent Hartree-Fock selfenergy - are well feasible [10] with the inclusion of correlations, simulations could 
only be carried out for very short times. This is a particular problem of ionization dynamics where one has to 
simultaneously treat bound and free electron states which requires large basis sets. 

On the other hand, the interaction of atoms with a strong laser field can, in many cases, be assumed to be a 
coherent process. Further, at sufficiently low density, the interaction between atoms can be neglected. Also, at 
room temperature conditions the density matrix of the atom is dominated by the contribution of the ground state 
wave function, and one can instead of quantum statistics apply a pure state description. 



2.2 Wave function based methods 

At the conditions described above the initial state \tpo) of an atom is given by the ground state wave function or a 
coherent superposition of atomic eigenfunctions. Then the dynamics following an external exciation is described 
by the time-dependent A^-electron Schrodinger equation (TDSE) with the hamiltonian (j2|i 

t^^\m) = Hit)\m), m = o)) = \^o). (6) 

The TDSE has been successfully solved in recent years by direct integration, exact diagonalization or various 
other techniques. For some recent results for atomic ionization with ultraviolet fields see [21-23] and refer- 
ences therein. Yet these kinds of exact caculations are presently limited to 3 ... 6 electrons, depending on the 
dimensionality. Therefore, approximate methods have been developed. A particularly successful approach is 
the Multiconfiguration time-dependent Hartree-Fock (MCTDHF) method due to Cederbaum [24], see also [25] 
and references therein. We have, therefore, implemented the MCTDHF method for photoionzation of atoms 
in a strong field and briefly outline the main idea. Starting point is the following ansatz for the many-electron 
wave function as linear combination of time-dependent Slater determinants which are denoted in terms of the 
occupation numbers rii , . . . um of the single-particle orbitals 

I^W) = H <^n(t) |ni,n2,--- ,nM]t). (7) 

n 

The sum runs over an integer number (^) of determinants. In the limit M oo the method yields the exact wave 
function, whereas for M = N one recovers the time-dependent Hartree-Fock approximation (TDHF). Thus, by 
increasing M one can systematically improve the treatment of electronic correlations. The equations of motion 
for the coefficients Cn and the single -particle orbitals \(j)a) are obtained by a variational principle and read, with 
a matrix notation (with respect to Slater determinants) 

iC{t) = H(t)C(t), (8) 



I 



\^°) 1 - I 0r )( 0r 



M \ M 



h \(t>a) + ^{P ^}^, Ptjlk Wjk \<Pl) 
ijkl 



(9) 



where pij and pij^i denote the single-particle and two-particle density matrix, and h is the single-particle hamil- 
tonian. Further, Wjk is the operator of the mean-field potential which, in coordinate representation, is given by 
VFjfe(r) — J df (pjiv) Eqs. ( |8|9| show that the expansion coefficients C obey a time-dependent 

Schrodinger equation, whereas the single-particle orbitals \(f>a) are the solutions of a more complicated non-linear 
integro-differential equation. 



3 IMCTDHF Results and Discussion 

We now present numerical results for photoionization of a four-electron one-dimensional model atom (beryllium, 
N = Z = 4) for which we presented nonequilibrium Greens functions results in ref. [10]. For the ID simulations 
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Fig. 1 Sketch of the ID 
Berylhum model atom, 
(a): Energies of the lowest 
Hartree-Fock orbitals and 
corresponding probabil- 
ity densities (vertically 
shifted), (b): total electron 
density prior to action of 
the laser pulse in Hartree- 
Fock (full blue line) and 
for M = 10 (dashed red 
line with higher peak). 



we use regularized Coulomb potentials, w(a;) — —Z{x^ + k^)^^^^, sind w{xi, X2) — [{xi — X2)^ + 
where the screening parameters are introduced to avoid divergencies. This model has been used successfully in 
many studies of atom-laser interaction where also the influence of k has been investigated. Following refs. [10,21] 
we use K = 1. 

A sketch of the atomic core potential and the (square of the) occupied HP basis functions (shifted by the orbital 
energies) is shown in Fig. [T] together with the total electron density. There are two bound states at the energies 
Eq = —1.371 and Ei = —0.312, each occupied by two electrons. The first task is to obtain accurate results 
for the ground state. The results of our MCTDHF calculations for the ground state energy are presented in the 
table. Convergence with respect to the correlation energy is reached for M = 9. In comparison, the ground 
state energy for the same system from NEGF calculations in second Born approximation is found to be equal to 
E2B ^ -6.7713, i.e. on the level of MCTDHF with M = 3. 



M 


2 


3 


4 


5 


6 


8 


10 


Energy [Hartree] 


-6.7395 


-6.7713 


-6.7800 


-6.7829 


-6.7847 


-6.7850 


-6.7851 



Table 1 Ground state energy of the ID Beryllium atom with n — 1 for different numbers of determinants M. We used a 
basis with 50 eigenfunctions of the single-particle hamiltonian which assured convergence with respect to the basis size. 



As a next step, we use the ground state wave function as the initial state and expose the atom to a time- 
dependent laser field of the form (|3]l with a fixed number of cycles choosing the pulse width according to t = 
107r/tj and an amplitude of £0 — O-l- We consider two frequencies, cji = |_Ei|/2andtj2 = (l-E'ol + l-E^i l)/2 which 
corresponds to a Keldysh-parameter 7 = Ip/2Up, [Ip is the ionization potential, and Up the ponderomotive 
potential. Up = £q/4uj'^] of 71 = 1.25 and 72 = 6.7. In the first case we expect that direct photoionization 
is not possible and the dominant process will be tunnel ionization. In contrast, in the second case multiphoton 
ionization should be expected. The time-dependent numerical results of our MCTDHF calculations are shown in 
Fig.|2] In the lower row we show the ionization yield [time-dependent expectation value of the number of ionized 
electrons Nion{t)]- Interestingly, Nion{t) is much larger for the smaller photon energy (case 1) which confirms 
our previous results [10]. This is due to the much larger electron excursion amplitude in the low frequency laser 
field [recall that, for a free electron, the amplitude is xq ~ w~^]. For the larger photon frequency uj2, xq is 
strongly reduced. Nevertheless, the two electrons in the upper state with energy Ei are driven substantially away 
from the nucleus during one half-cycle of the field but they are driven back during the following half-cycle. As 
a consequence, the ionization probability is substantially reduced. In this second case, MCTDHF calculations 
with different numbers of determinants M are very close to each other, and even a TDHF calculation yields the 
qualitatively correct result. More interesting is the behavior for the smaller frequency. Here two electrons are 
driven far away from the nucleus and they occupy continuum states which leads to strong correlation effects. As 
a result, TDHF calculations - while overall giving the correct picture - become inaccurate already after about 
one laser cycle (t > 40). 
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Fig. 2 (Color online) 
Ionization dynamics of a 
ID Beryllium model atom 
for two laser frequencies 
and different approxima- 
tions for the electronic 
correlations. Left column: 
u] = 0.156 a. u., Right 
column: uj = 0.842 a.u. 
Top row shows the 
mean electron coordinate 
and bottom row the 
mean number of ionized 
electrons. 



Summarizing, we have presented MCTDHF calculations for a four electron model atom exposed to strong 
laser pulses with two different frequencies. These calculations are rather efficient and allow to follow the time- 
resolved ionization dynamics and the history of electronic correlations in single atoms. They show that, overall, 
TDHF gives the correct tendencies. However, it fails to accurately describe the dynamics of electrons at longer 
times, in particular in cases of strong ionization. In future work we will concentrate on correlation dynamics and 
on the inneratomic relaxation as well as on further comparisons with Nonequilibrium Greens functions results. 
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